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seasonal tl'.ernocline is acconpanied by iiany rand on fluctuations. The 
physical processes causing the decay and the fluctuations in the nixed- 
layer depth are cxai.uned and qualitatively evaluated, 

Tne Tukey spectrun analysis progranned for the CDG 1604 electronic 
digital computer was used in analysis of oscillations at the bottom of 
the mixed layer to determine the distribution of \;ave energy for series 
of data taken at hourly intervals. -An energy peak centered neai- 12 hours 
was obsei’\’’od to prcdoi.ainate, and this energy is equiva.lent to that of a 
12-hour internal \;ave \;ith height of 27,0 foot. The computer \;as also 
used to select the significant meteorological and oceanographic parameters 
which could be used to predict a change in the mixed “layer depth by iise 
of the BUD 07 multiple regression proqprrm. The most significant paa’ameter 
was sea surface temperature. The best correlation coefficient for this 
parameter occurred for lags of aero to 12 hours. 

Uie results of the analysi s support the theory that convection is 
the physical process which causes the seasonal decay of the thermo- 
cline dtiring the cooling season and that short-term fluctuations of the 
mixed-layer depth are due primarily to internal \;aves. 

The authors wish to orpress their appreciation for the assistance 
and encouragement given them by Associate Professors G, II. Jung and 
J. D. Wickham, and by Professor D. A. Williams and the staff of the com- 
puter center at the U. S. Naval Postgraduate School in this investigation. 
The authors also wish to thanl: the Pacific Oceanographic Group of Canada 
for malcing the oceanO;p?aphic data available. 
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1 . 

The vertical tenporature structure of the near- surf ace layer in the 
open ocean is a conplicated phenomenon which varies in space and tine. 

Certain aspects of the temperature structure arc of great significance 
in determining sonar effectiveness. Thus, both the military and the com- 
mercial fishing industries are interested in prediction of the thermal 
structure. 

Various processes continually modify the thermal structure to produce 
an essentially homogeneous (mixed) surface layer, underlain by a thin 
layer of rapidly increasing density, v;hich varies in thickness through- 
o'af the year, Tlic aixed«»layer depth (referred to as the I-ILD) is the depth 
belov/ the v;atcr surface to which mixing has established isothermal condi- 
tions, The lov/er boundary of the IILD is the top of the thormoclino, a 
thin layer of largo, negative vortical temperature gradient, usually associa- 
ted with the layer of increasing density, A complex interaction of meteoro- 
logical, oceanographic, and (indirectly) astronomical factors provide the 
driving forces for changes in the ocean thermal structure. 

Random and periodic fluctuations of the KLD have been examined by 
many investigators. It is knovm, as a result, that the iILD may vary from 
a fcv7 inches in the heating season to several hundred feet in the cooling 
season. At times it may fluctuate tens of feet in a matter of minutes. 

The objective of this research is to identify more clearly significant 
processes which affect the IILD during the cooling season and to relate the 
IILD to meteorological and oceanographic parameters. The processes capable 
of affecting the vertical temperature distribution are discussed. in the 
appendix. Certain of these processes will be considered in this paper: 

(l) thermohaline convection, (2) dynamic convection due to wind and vavo 
action, and (3) unidentified processes which generate internal waves. 
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Up to now, unly yaalit t_vj uoocri^jtiono of the ..axod-loycr behavior 
durin{:; the cooling season are foiind in the litcratui-e. These usually 
state that increased wind stirring and surface cooling act together during 
the cooling season to drive the seasonal thcrraoclinc do\.'n to its naxinun 
annual depth. Our findings will shovr that, at Station Paoa, thernohaline 
convection is the nest important factor and that v;ind nixing (dynamic con- 
vection) is actually insignificant during the cooling season. 

Met e Orel ogical and oceanographic pararacters contributing significantly 
to the change in the MLD during the cooling season arc deterninedj and from 
those the physical processes of importance arc inferred. The investiga- 
tion was carried out with the aid of statistical correlation and wavo- 
analysis programs, the BIID 07 [l^ and the lUkcy spectrum analysis [ 153 , 
utilizing the GDC I 6 O 4 digital computer. Regression equations and correla- 
tion coefficients arc obtained by the BIID 07 program. One part of the 
Tukey spectrum analysis was used to determine cross spectra of the IDD 
fluctuations and the sea-surface temperatures. This program was also used 

t 

in the analysis of internal v;avcs. 

The investigation is ba.scd on oceanographic and weather da.ta obtained 
by the Fisheries Research Board of Canada, Pacific Occanogi’aphic Group, 
for the months of October, flovcmber, and December during the years 1956, 
1957 and 1958 [6. 7 ]. 
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2 . '3ackf,Tound 

The attempt to identify tlic factors v/hich bring about the ho\irly 
fluetuations in the lILD v;as based on the authors' earlier endeavor to 
derive an empirical method to predict the average HI.D on a daily basis 
during the cooling season. It became apparent then that the hourly fluc- 
tuations of the I ID could be of equal or larger magnitude than the average 
deepening of tlio tlLn during one or more days; thus, to use one or an 
average of tvra or more observed I>XD values per day as representative of 
the depth for any particular hour of that day is misleading Q J . The 
probable cause of these fluctuations is internal v;avcs. 

Students at the U. S. Naval Postgraduate School during 1961 and 1962 
derived empirical methods to predict the M'D at Ocean Station Papa dui’ing 
each of the 12 months of the year. In addition, the onset of the seasonal 
thcrmoclinc at Ocean Station Papa, uhicli talces place in April or Kay, v:as 
studied by Clark in 1961 [2j* Clark c.uccessfull.y related the thermocline 
onset to upper-air parameters. During the heating season, v;hen the sea- 
sonal thormoclino gradually increases in magnitude, vjind mixing processes 
have been clearly established as the dominant factor in dotoiraining the 
I JiD, by Tabata [22^ end Geary [lo]] in 1961. The cooling season remains to 
be studied. V/e have made this attempt. 



3. Dcccrijition of Ocoan St-ition Papa 

The geographic location of Ocean Station Papa is 50N, 145W (figairo l). 
It is approximately 400 miles north of the boundary botveon the eastern 
subarctic Xv'ater roass and the subtropic vator mss [8] . The vrater depth 

at Papa's location is 2000 fathoms; therefore, the therml structure is 
in no v/ay influenced by bottom topography. Also, its location is approxi- 
mately in the center of the A].askan Current G;/ral Q 9 ] (figure 1 ). 
Because of this location, the effects of advcction arc minimized apd will 
not be considered. 

The weather ship normally maintains an on-station position \;ithin 
a ter— mile square centered on 50H, 145W, In order to avoid rough weather, 
the ship may steam anyv/hore \;ithin the 210-milc square grid (figure 2), 
Sach bathythermograph (3T) sounding is identified by a tvro-lotter code. 

If the ship is on-station, the letter group OS indicates the position. 
The ship position for the greater part of the period of our concern v;as 
v/ithin a 30-r,ilc square centered at OS. .Iny data outside of this 30-milo 
square were not considered in this study. 
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oi To;..' •'r •. 1..-’. S'llinity oiructu:”c dvi’inj tl'.o Cooling Sc.inon 
x'^iguvo 3 ro,,rooontr; tho docny of the th3rr;oclino during the yoor 
1956 . In A’agunt, the rurfacc Inyoro have boon heated to a mxinura, approxi- 
mately 37F. Due prinai'ily to wind mixing, there is a sliallou nixed layer 
approxirately ICO foot in depth with a sharp underlying thernoclinc. With- 
in the thcrnocline, the water temperature drops to 4 OF at a depth of 200 
feet. Progressing into the fall and vdnter, the M-D continually increases 
and the magnitude of the thernocline docrcasoc. By the end of December, 
tlae MTjD is nearly 330 feet in dcptl:^ and the sea-surface tonpcratui’o has 
uccroasod to approximately 42F. Due to the presence of an intense halo- 
cline, the MIjD I'cach.os a limiting depth of 330 feet in Jantiarj^ and iso- 
thermal conditions' prevail down to this depth for the remainder of the 
winter. 

The basic salinity striicturc for the eastern subarctic Pacific v/atcr 
mass is shown in figure 4. Tliroc zones arc indicated. The upper zone 
extends from the sea surface to a depth of 320 foot, and is characterized 
by rolatively lo\;-salinity v;ator (32, '59^ f2l] la the fall, the \;ppcr 
zone is isoh.aline and is narlceu by the presence of the seasonal thermo- 
clinc. Dccaixso of the Isohallnc conditions in the upper 330 feet, vre 
assume that significant variation in density is governed by changes in 
temperature only. 

The haloclinc represonts a transition zono botv/oon the upper and 
lower zones. Kero, the salinity increases markedly with depth (as much 
as one part per thousand v.’ithin .on interval of 330 feet). The haloclinc 



represents a limiting depth for the seasonal thormoclinc. It is similar 
to having an ocean bottom at this depth as far as convective nixing is 
concerned. Tlxe density increase is so great through this layer that mic'ing 

cannot take place through it ^20] . Within the lower zone 0^50 feet) 
the salin i, -"adually incroac 



cs xdth depth to the ocean bottom, 
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Tho purpose of this statistical study of the MLD is three-fold. 

All attnnpt is mde: (l) to dctorirdne the relative inportanec of thcrr.o- 

halinc conveetive mixing and wind-inducod mixing on tho deepening of tho 
thormoclino during the cooling season; (2/ to determine tho lag betv;oen 
a signifioant meteorological or oceanographic event and a subsequent cliango 
in the MLD; and (3) to examine tho inportanoo of the seemingly random 
short-term fluctuations of the MTjD and to provide a description of those 
internal v;avcs present. 

First, the parameter or parameters are identified which have the 
largest linear correlation viith the deepening of the MiD. Next, a further 
investigation has been undertalcen of these significant parameters to deter- 
mine the lag in the mixing process, that is, the time between the onset 
of a change in the parameter and a subsequent deepening of the MLD. 

Finally, the observed variations in the M.D v/hich eannot bo aeoounted for 
by either mixing process arc investigated, 

5,b, Method of Investigation 

The parameters to bo investigated are sub-divided into three groups. 
Those related to convective stirring arc: (l) temperature of the sea, 

at tho time of the sounding and v;ith a 2-hour lag; (2) air temperature, 
at the time of tho sounding and v;ith a 2-hour lag; (3) devupoint tempera- 
ture and v;ct-bulb temperature, at the time of the sounding; (4) relative 
humidity, at 2-hour and 5-hour lags; and (5) low-eloud cover, with a 2- 
hour lag. 

Parameters chosen to represent v;ind and v/ave mixing are: (1) v;ind 

speed, at 2-hour and 5-hour lags; and (2) height of significant waves, 
at 2-hour and 5-hour lags. 



Parameters chosen ’.,’hich verc chov;n to have a correlation v;ith the 
th.crnal structure by others in previous studies, but did not conveniently 
fall into the abo'^’c categories, are; (l) 500-mb height, in goopotential 
meters; (2) 500-rab temperature and 500-nb relative hixmidity, at a 2-hour 
lag; and (3) surface-pressiu’o change, v;ith a 2-hour and 5-hour lag. 

Those data were obtained from synoptic 3-hour ly reports, tv,’ice-daily 
constant— pres sure records, and the BT traces. The available BT‘s for the 
period of investigation wore tabulated in tvjo groups, 0200S and 17003, 
for each year to eliminate any possible yearly anomalies or diurnal trends. 

An investigation by Soripps Institution of Oceanography of the annual MLD 
fluctuation in the northeastern Pacific Ocean shows that the MLD most jfre- 
quontly observed in a given month often varies oonsidorably from year to 
year P8]. The available literature on diurnal trends during the cooling 
season is not so clear. A preliminary investigation of the data showed 
no marked diiornal trend. This indication v/as supported by Dr. Tully of 
the Pacific Oceanographic Group, in a conversation with the authors, in 
v/hich he stated thau during the cooling season the MLD is too deep to have 
a clearly defined diurnal fluctuation, 

With the aid of a multiple-regression and correlation analysis, the 
BIMD 07 programmed for the GDC 1604 high-speed eomputer, the linear- 
correlation coefficients were obtained for the prc-solcctcd group of me- 
teorological and oceanographic parameters with the MLD. The BIMD 07 pro- 
gram was used because it can select different sub-samples of data that 
are obtained from the sane population. Thus, the program would perform 
linear-coirolation analysis and linear mltiple regression on the data 
for 02003 and 17003, separately and in combination. In this v;ay, a close 
examnation of possible diurnal trends could be made for each of the throe 
years studied, 1956 through 1958. The analysis provided regression equations 
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for nine difforent data groups for the periods shown in Table I. 

The BII'D 07 program yields best results v;hcn the nui.-iber of indepen- 
dent variables is not greater than one-half the number of values of the 
dependent variable. Since the smallest number of values for the depen- 
dent variable is 36 (at 02003, 1958) only 18 independent variables were 
investigated. A set of 18 secondary parameters \;as tabulated in the event 
that the analysis of the 18 primary parameters proved inconclusive. 

Tlicrc was no attempt by the authors to eliminate possible duplication 
of parameters that arc indicative of the same process. For example, dc\-^- 
point temperature, wet-bulb temperattirc, relative humidity, and low-cloud 
cover are all related to evaporationj but it was left to the BIHD 07 program 
to discern which -ones, if any, were significant in contributing to the 
prediction of the deepening of the liLD. 

The BIM) 07 is unable to produce a multiple linear-regression equation 
if any independent variable is a linear combination of tv/o or more inde- 
pendent variables. Therefore, temperature of the air minus temperature 
of the sea, \jhich is often used as an indicator of thcrmohalinc convection, 
could not be investigated in conjunction vzith the air temperature and sea 
temperature. Consequently, three additional problems were computed using 
the air-minus-scc fcenperature parameter in place of sea-temperature para- 
meters. See Table II for the tabulated results of the 12 problems. 

Tn analysis of all the data, certain facts and tendencies ore evident; 
these arc discussed in the following two paragraphs. 

The only parameters that shov; a consistently high linear correlation 
with the MT.D in every sample or sub-sample arc the sea-surface tempera- 
tures, at the time of g sounding and v/ith a 2-hour lag. They also accounted 

"*(Tliis small num.ber results because BT data for October, 195§ arc miss- 
ing; only a small nmber of observations was available in 1957 as vjcll, since 
December, 195*^ observations were not taken. ) 
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for alnoGo all of tiio varianco explained by the regreGnlon equation using 
all 18 paraiaetcrs. Figure 4 shows the close relationship bctv/cen the 
sea- surface tenporature and the MLD from October through December, 1958, 
Likcvfisc, Zcttcl [ 25 ] suggested that sea- surface tomnerature might be a 
good "indicator" for the IXD. 

V/ind speed and height of the significant waves, the parameters chosen 
as indicators of wind stirring, have either very lovr negative or lov; posi- 
tive correlation coeffieients. 

In order to determine the time required for an increase in convoetivo 
activity due to a reduction of sea-surface temperature to effect a sub- 
sequent deepening of the MjD, the llLD’s for 17002, 195 were tabulated 
v;ith the sea temperature at lags of zero, two, five, eight, eleven, four- 
teen, seventeen, tv/enty-onc, and tv;enty-three hours as the independent 
variables. A simple linear-eorrolation prograra was used to determine 
that the lag yielding the best correlation was two hours (Table III). 

Also, two one-hourly series of KLD for 79 hours in October, 1957, and 
for 51 hours in November, 1958, v;cre tabulated v;ith their eorresponding 
sea-surface temperatures. The energy-density spectrum and the eross-spcc- 
tnim of both the IITjD and the sea- surface -temperature time series i:as ob- 
tained from the Tukey analysis. Maximum lags of 39 hours and 24 hours 
v;ere used for the 79-hour and 51-hour series respectively. The fv/o co- 
spectio for the 79-hour and the 51-hour series show largo values for 
h’s around 8 hours and 4 hours respectively (figure 5)» 

From the results of the simple linear correlation and the Tukey analysis,^ 
it is soon that all lags from two to six hours show high correlation. 8c- 
yond 12 hours the correlation rapidly decroaser, 

A linear-regression equation was derived for 17002, 1956, using sea 
temperature v;ith a 2-hour lag as the single independent variable (Table IV ). 



1 on cross-correlation, not shown 
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Tho fact that the results of this equation as measured by the standard 
error conuarc favorably with tho multiple- regression equation derived from 
13 independent variables for the same year supports the initial analysis 
of the problem. If a regression equation vjsrc to be used to forecast 1I.D, 
one using only the sea- surface temperature taken with a lag of up to 12 
hours would be much more reliable than an equation with more parameters. 
This is true simply because it would have a larger number of degrees of 
freedom. Note that all of the regression equations have a standard error 
of from 15 to 30 feet, with an average for all problems of 22.15 feet 
(Table V, Table VI). 

There are several possible explanations for this standard error of 
more than 20 feet: (l) ft is possible that not enough parameters v/crc 

used, and that by substituting for primary parameters v/hich sho\/ed a very 
poor correlation, some of tho secondary parameters or primary parameters 
with different lags, the standard error of the new regression equation 
would be reduced; or (2) there is a strong possibility that a non-linear 
relationship between a parameter and the MLD could account for a large 
part of the standard error (hov/cver, the BIMD 07 program cannot detect 
this); or (3) the error may be due to internal waves. Of these possi- 
bilities, (3) seems tho most likely reason to explain tho observed stan- 
dard error. 
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6. ilnalysis of Intornal VJavoc 

The base of the jaLxed layer in the ocean is cubjcct to nany vertical 
oscillations. Figure 7 sho\Js a series of daily MLD's for the period Octo- 
ber through Dcccnbor, 1956. Each plotted point on the curve is an average 
daily value of the IXD. Tv;o individual soundings arc normally used; how- 
ever, up to 24 soundings wore available on several days. It is observed 
that the MIjD undergoes a seasonal downviard trend. During the October 
through December period, the seasonal trend of the MjD finds it increas- 
ing from approxiriately 130 feet to 330 feet. Tliis is an average deepen- 
ing of 2.17 feet per day. Also, largo vertical fluctuations are super- 
imposed on this dovmuard trend, v/hich become larg«‘" in magnitude as the 
nixed -layer thickness increases. Note in figure 7 that the vertical fluc- 
tixations of the ILD can be ten times greater than the average daily deepen- 
ing. In addition, it is seen (figure 8) that hourly fluctuations can bo 
of equal or greater magnitude than the average daily deepening during the 
fall. 

Figure 8 represents a one-hourly tine series for 51 consecutive hours 
during November, 1958 (series l). A second one-hourly time series rep- 
resents 79 consecutive observations taken during October, 1957 (series 2). 
The magnitude of the IILD fluctuations is about the same, whether it ir 
computed from hourly, twice-daily, or daily observations. The average 
fluctuation per hour of the KLD for series 1 is 18 feet. The average 
fluctuation between the consecutive 02003 and 17003 soundings for the month 
of November, 1958^ is 21 feet. Also, the averaged daily fluctuations during 
November is 25 feet. The nurabers for series 2 arc similar. A December, 
1956 , one-hourly tine series (figure 9) Was available but was too short 
for detailed analysis. Hov;ever, it illustrates the extreme fluctuation 
of the I'iLD v;hich talces place over a short time interval. Over a 3-hour 
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,)criod, tl. . 1 LD incr'’'ircd fror.i 270 To t to 370 foot, and then decreased 
attain to 235 foot. 

Since nixing processes can account for only an increase in the IXD, 
and since observed fluctuations are in both directions (up and do\;n), 
and arc a consistent feature of the data, these observed fluctuations 
arc considered to bo due to internal vaves. 

Internal ;;aves nay occur \;ithin stratified water and in water in 
\;hich the density increases \iith depth. The largest vertical displace- 
nents of tlie \ratcr particles arc to be found v/ithin the boundary between 
layers of different density. The displacement ai:iplitudo diminishes above 
and bclov/ this boundary, approaching zero at the sea surface and the ocean 
bottom. The density difference betv/een the waters above and below the 
seasonal thcrmocline is, of coiuse, much less than the density difference 
between the air and the v;ater at the sea surface. Therefore, a boundary 
v/ithin the ocean can be much more easily displaced than the surface of 
the sea CL7] Ufford [24] states that surface v/aves require 30,000 times 
more energy to start and maintain than is required for internal v/aves of 
the same amplitude. Energy for internal v/aves could come from surface 
disturbances or from currents v/ithin the ocean. 

Besides rapid vertical oscillations in the IIjD, it v;as observed that 
the thickness (or vertical extent) of the thcrmocline layer bclov/ the 
lILD may vary markedly v/ith time. Figure 10 shov/s the thcrmocline thick- 
ness for scries 1 . Variations in the thcrmocline thickness indicate that 
intcrnal-v/ave characteristics of phase and amplitude vary v/ith depth, 

Mso, the temperature gradient in the thcrmocline varies markedly because 
of changes in its thickness. The actual temporature diffcrcnc.: tlirough 
the thcrmocline is a relatively stable feature and is uh-'rafore predic- 
table jll] . 
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loy n, ;.~rl 'i ■,.. .. V.-? ^ ir.Vctjc co\’J-d fcr). :ind ,;roya3'',to di’rin,'; t!:.:; nui.mor 



at the level of tlij eeaconal ther:.iocllno at Ocean Station P^pa. Vj'ith 
one-hourly oboervationo available for analysis, the :.iininun neried that 
could be Invortipated is t\:o hoairs. 



The 'Pukey spcctrun analysis, proyraairied for 
i.'as used to analyse each tine series. The min 
eredoriinant periods and •iscoclated anplitudes pr 



the 3 DC lb04 ccr.iputcr, 
object y'as to find any 
esent in tl:c IXD oscil- 



lations, The Tukey analysis cctimtos the enor^y-density as 
of i^oriod in a yiven series. This is done by conpi'.tiny tlio 
tion function for a nur.ber of yiven lags in t!'.o series; then 



a function 
autc-corrcla- 



a snoothed 



spectral ostinato is obtain-^d 



fron the Fourier ccsinc transform, of the 



auto-correlation function. 

Initially, a correloprar ’..'as plotted for both time scries (figure 11), 
The corrclocra::i is a plot of the non-nornaliacd a’ato-corrclation coeffi- 
cients as a function of the lag (ic) , Inspection of this corrologr<ar.i appears 
to indicate th.at th.c spectrum of series 1 has a relatively Large energy- 
density mrcinuji near the 12-hour period and a secondary maxii.juri vrith a 
peak near six hoars; and that the spcctnnn of series 2 has only one naxinun 
near 12 hours. 



Kext, the energy-density spectrum of both tine scries v?as obtained 
from the Tulcey analysis. Figure 12 i s a plot of the energy-density spectrum 
for series 1. A maxinum lag (n) of 24 hoixrs \;as used, uhich alloyed the 
spectrum to include a range of periods from tvo hours through 4S hours. 

The spcctryn helps confirm, the indications of the corrclogram. A ma:cimura 
of energy of the variance of tho liT.D is sho\m to be produced by fluctua- 
tions with periods on the order of 12 hours. Also, a large energy-density 
is indicated for periods around six hours. There is an energy gap for the 
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periods eround eight hours, and the energy for periods less than five hours 
drops off rapidly to oji insignificant ai.iount. 

The total energy is a known quantity, equal to twice the variance. 

By measuring total area under the spectrum curve and the area under each 
frequency peak bctvjcon half-energy points, equivalent amplitudes for v;aves 

a 

\7ith periods near 12 hours and six hours vroro computed. These computed 
amplitudes arc 18 feet and 12 feet for the 12-hour and the 6-hour periods, 
respectively. 

The energy-density spectrim for scries 2 is shovm in figure 13« A 
narciraun lag (n) of 30 hours \-7as used, which allovied inspection of periods 
ranging from tv;o hours through 60 hours, A relatively large energy maxi- 
mum is again indicated for periods around 12 hours. The curve then falls 
off rapidly, but shov;s a slight maximum near the 6-hour range of periods. 

A single sine v:avc with energy equivalent to that under the 12-hour pealc 
would have an amplitude of nine feet. 

This analysis reveals three features of the thermal stiaicttu'c at 
Ocean Station Pana during the tine intcrva.1 investigated, 

(l) Vertical oscillations of the MLD are present at the thcrmoclinc 
during the cooling season. Tne magnitudes of these oscillations are 
similar in all data, \;hcther they be from hourly, txiicc-daily, or daily 
obsor'/ations. 

^2) Although the time series available \;as much too short on which 
to base strong conclusions, the analysis did reveal that the raajor energy 
peaks occurred \Jith a period of approximately 12 hours vrith an average 
amplitude of 13.5 feet. 

( 3 ) The value of the standard error of estimate resulting from the 

BUD 07 regression equations is very close to the height of the postulated 
12-hour wave (that is, 27 feet). Therefore, it is concluded that there 
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ic a direct relation botuoon this error and the presence of internal 



\:avcs during; the period of investigation. 
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7, Conclunionc and Rocoiicacndationn 

As a result of this study of the influences on the TiLD during the 
cooling season, the follo\;ing conclusions can be reachea. 

(1) Thernohalinc convection is the important process for the deepen- 
ing of the ML'D in the autumn season. 

( 2 ) Dei ping the period of study, the M.D is too deep to bo affected 
by v;avc stirring or diurnal trends. 

( 3 ) Th2 tine lag bct\;cen surface cooling and deepening of the ILD 

by the convective process is snail; this caiggests that once the cooling 
season begins, convective circulation is essentially continuous. This 
circulation nay be maintained during periods of near-zero heat across 

the air-sea interface by wind energy, as suggested by I.aFond [«] ; hov:- 
cver, \mtil an injection of potential energy in the form of rapid surface 
cooling, the M.D vill remain relatively static. If the increased surface 
density caused by sui'face cooling is sufficient, the convective circula- 
tion intensifies and drives the MLD deeper. 

( 4 ) A major portion of the MjD variability appears to be associated 
■v.'ith internal \/aves. The analysis of internal v/aves revealed that the 
raajor energy pealcs occurred with a period of approcinatcly 12 hours v;ith 
an average amplitude of 13.5 feet, 

Karly in 1953, at a Navy- sponsored thornoclinc conference attended 
by 20 of this country’s most prominent oceanographers and meteorologists, 
the follo\.’ing question uas asked his colleagues by Dr. T. F. Malone; "How 
accurate a prediction of the MLD and its hourly fluctxiations is required 
by the military?" [ 19 ] This remains a vital question, 

Tf the military is \/illing to accept forecasts of the MLD within 
ranges of 20 or 30 feet, then existing methods of prediction will suffice. 
Any increase in accuracy over existing schemes for predicting the MLD at 



16 



a n'’,i’ticul’r lioiu' on a narticular day can bo acconpliohcd in aroac vrhorc 
internal viavcc aro inportant only after prediction of tho internal \:avc 
pattern is nossible. This is a conplicatcd problen. 

If a sufficiently lon^, detailed tir.ic series of BT's could bo mdc 
available, intornal v;avcs could be filtered out by use of a technique such 
as proposed by Linnettc [u]. This uould leave a tine series of the M.D 
that is due entirely to convective miring or other causes, thus enabling 
a mch better quantitative analysis of tho convective process. This is 
suggested as the next logical extension of the research reported in this 
paper. 
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Figure PoGition-Indicating Grid for Ocean Weather 
Station Papa, vith a Ilercator Projection of 
a Lao-tude and Longitude Grid Superinposod 
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Fi^xiro Groc 33 “-Spec - XD and Sova-Svirfaco - 

Tonporaturo Tine Series (llovoraber^ 19^S) 
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Figure 7. Average Daily MLD 

(October through December, 1956) 
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Figure 11. Correlogram for One-Hourly Time Series 
(October, 1957, and Novembfer , 1958) 
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Table I 





1956 


1957 


0200S 


43 




17003 


61 


46 




109 


90 






(02003 i 
(17003 i 



19^8 

36 

— 

79 

1600 local noan tlKo) 
0700 local ncan tine) 



Bathythornograph Soundings Available 
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I 




T:'J)1e II 







02005 






Par-uietcr 


1254. 


1957 


_li'48 


I9i6 




-.905 


-.903 


-.649 


-.850 


Ta 


-.635 


-.583 


-.354 


-.410 


Ta-Tw 


.369 






.372 


^dp 




-.146 


.104 


-.162 


Tv/b 


-.519 


-.284 


-.008 


.153 


Ril 


.167 


.266 


.272 


.121 




.028 


.266 


.333 


.069 


F 


.180 


.107 


.212 


.351 


FF 


-.046 


.339 


.012 


.110 


'V 3 


-.200 


.686 


-.089 


-.077 


(Tw)-2 


-.925 


-.879 


-.793 


-.884 


(Ta)-2 


- . 6S4 


- . 646 


-.228 


-.382 


(Ta-Tw)-2 


.314 






.448 


(FF)_5 


.019 


.302 


-.073 


-.262 


(Hi/3)-5 


-.200 


. 6S3 


-.037 


-.223 


P 


-.311 


-.134 


.125 


.039 


(P )^5 


-.060 


-.004 


-.101 


.041 


GPM5 


.118 


-.463 


-.043 


.215 


'^5 


.113 


-.439 


.091 


.166 


RIT5 


-.136 


.285 


.150 


.137 



1700S 




020C5 and 170CS 


1957 


1418 


1956 


1417 . 


1958 


-.767 


-.814 


-.871 


-.829 


-.734 


-.386 


-.331 


-.502 


-.489 


-.333 






.369 






-.020 


.017 


-.242 


-.050 


.062 


-.189 


-.115 


.092 


-.211 


.063 


.290 


.202 


.1/a 


.263 


.238 


.400 


-.015 


.049 


.327 


.148 


-,.093 


.249 


.276 


.004 


.231 


.307 


-.131 


.047 


.360 


-.062 


.487 


-.058 


-.127 


.581 


-.060 


-.649 


-.843 


-.899 


-.726 


-.819 


1 

1 


-.314 


-.505 


-.502 


-.277 






.383 






.352 


0,172 


-.143 


.311 


-.115 


.525 


-.389 


-.192 


.583 


-.245 


.170 


-.342 


-.107 


-.027 


-.090 


.327 


-.215 


.047 


.125 


-.158 


-.419 


.143 


.169 


-.437 


-.050 


.623 


.080 


.U5 


.181 


.080 


.045 


.186 


.028 


.166 


.154 



Sinplo Lincar-Gorrclation Goofficicntc 
botv/cen the IS Primr .7 P-’r-inotcrc! and the HD 
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Table III 



Soa-S’oi’f ac o- reMpr-ratui’G 
Lar in Koiirs 

0 

2 

5 

0 

o 

11 

1/t 

17 

20 

23 



Coi-rolatioa Coofficicn 
w itli the IliD 

-.04975 

-.88409 

-.87828 

-.37576 

-.37173 

-.86479 

-.86829 

-.86125 

-.85189 



Sinplo Linoar-Corrclation Coefficients 
for Sor^Surf-'LCo Te. loeratures and the IXD 
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Table IV 



Linear-correlation coofficiont 
(for To_ 2 ) -.S 84 

Standard cri*or of the estimte 30.33 

Range of residuals 121.2 

Regression equation; = 1152.94 - 20.24 (Tj,_ 2 ) 



Regression Results for One Independent Variable 
(17003, 1956 ) 
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Table V 



1956 


0200S 

1957 


1958 


1956 


1700S 

1957 


1958 


02002 and 17002 
1956 1957 1958 


(Ty and 


T ^,;_2 included, T 


a-Ty and (Tg^-Ty). 


.2 onitted) 




.950 


.941 


.864 


.915 


.831 


.915 


.909 


.874 


.863 


26.5 


14.4 


26.7 


31.4 


18.0 


19.4 


29.9 


15.8 


21.0 


80.4 


45.4 


29.9 


104.8 


49.6 


59.8 


118.8 


57.1 


102.5 


210.1 


179.4 


273.3 


212.7 


170.5 


279.6 


211.6 


174.8 


276.8 


67.3 


32.4 


37.0 


64.9 


25.0 


36.27 


65.7 


29.0 


36.5 


and (T^- 


•lvr)_2 included, 


Ty and Ty. 


onitted) 




.950 






.917 






.910 






26.5 






30.9 






29.8 






80.4 






104.4 






121.8 






210.1 






212.7 






211.6 






67.3 






64.9 






65.7 







Ihltiple-Regrossion Recults for 18 Independent Variables 



Table VI 



I'iultiplc-sorrelation coefficient 



(for 18 independent variables) .896 

Linear-correlation coeffici ''nt 

(for Tg and Tg_ 2 ) -.814 

Standard error of the estimte 22.15 

Ranf^e of residuals 72.0 

Standard deviation 43. B 



Thi'oe-Year Avoi'af^es for Regression Results 
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Heat T ransnort in t’no Sea OJ^d ?.xc hanf :e acronr. ^\q Air- Soa In terface 

A. Thernohaline convection . The vertical displacencnt of orall oar- 
cels of vater occurs U’hon a srall part of n v;ator nass is heavier than the 
i,;ator underneath it. To restore this unstable condition to eq\iilibrivm, 
the heavier water tends to sink \/hile the lighter v/ater rises. Associated 
\/ith these forced vertical novements of snail parcels there is also a trans- 
port of the chara.cteristic properties of oca water in a vortical direction 
which leads to an equalisation of any vertical differences in these propci*- 
tlco 'Jhich my be present. It is believed that convective ni:cing is as 
iraportant an agent of heat transfer in the ocean as wind stirring during 
the cooling season ^2, 12. 2(^ , 

An initial increase in the density of smll particles at the s^orface 
accorapanics an increase in salinity (due to evaporation, or to the form- 
tion of ice) or a decrease in temperature. One or more of these in combi- 
nation my be involved in thernohaline convection. In lower latitudes, 
v/here there are only smll variations in the temperature, the heat loss 
is outweighed by the effect of evaporation; in polar regions, in addition 
to radiation and evaporation the increase in salinity due to formtion of 
ice is also effective. In temperate latitudes the heat loss by radiation 
is the decisive factor. 

Convective sensible heat exchange (conduction) between air and sea 
surface has been studied by Kuhbrodt [A] for the North Atlantic. He has 
shoim the convectional heat flux to be anproximatcly 20 gm cal/cn^ per day 
during the period of the year when the air is cooler than the sea surface. 
This value, oven as a rough estimte, is too large to be neglected in the 
heat budget of the ocean. If the surface of the sea is warmer than the 
atmosphere, the air is heated at the interface and the vcrt3c.-xl stratifi- 
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cation of the air becomes unstable; the air bccor.es turbulent -./ith a con- 
stant replacement by cool air at tire sea surface and the vertical heat 
transport becomes large. 

Values given by Defant for the heat loss in \inits of fjr ral/cm^ 
per day for 50N latitude are 116 units duo to effective back radiation, 

78 units duo to evaporation, and 20 ixnits due to convection of sensible 
heat (conduction) as a yearly average. Certainly, during the cooling 
season heat loss due to conduction vould be even more significant. These 
values are derived, assuming that the heat exchange through the ocean sur- 
face occurs independently for each separate latitude belt. Therefore, no 
meridional heat exchange (by ocean currents and by horizontal nixing) \;as 
considered, 

B. Dynami c convection ( forced vertical nixing due to w ind and wave 
action ) . This process has been v:ell studied and considered to bo the 
dominant process of vertical nixing in the heating season. Laevastu (l?^ 
has proposed an empirical formula for determining the MLD by significant 
v;ave height alone. 

C. Dynamic convection due to forced vortical nixing of ocean currents 
can result in convective nixing, and influences the ME, Tliis 
process \/as not considered because of the difficulty of measurement; hov/- 
evor, since the Station is located in the center of the Alaskan Gyral, 
ocean— current effects are probably smll, 

D. Therma l conductivity can take place ’..'hen a ^'ortical temperature 
gradient exists in the ocea:' , Heat is transferred by molecular heat con- 
duction processes. Defant |4j concludes that this process is insignificant 
in oce.anogranhic investigations due to the extremely long period involved 
in the conduction process. Therefore, this process v;ill not be considered 
in this study. 
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^ ntrai!i:nont d~'o to turbulencs at tl’o Ycnoclinc ♦ Durin;;; the 
coolin;' acason in the open ocean, a v/ell-clefined ioothornal surface layer 
exists v.’hich e::tcnds dovm.'ard to an interface of rapidly decreasinp; ten- 
)3oratui'c, Investi';;ationG of the subarctic water indicate an isohaline 
condition to a depth of 331 feet, so the effect of salinity on density 
can be ignored [ 23 J , and thus \io my consider the therriodine as a pycno- 
cline. Because less dense waruer v/ater overlies dense cold \:ater, v/e have 
a stable condition mch like a surface inversion in the atr.iospliere, When 
surface cooling takes place, thereby crea.ting an unstable condition, free 
convection can occui’. Cronwcll|”3j has denonstrated that this free convec- 
tion in a nodel causes a turbulent exchange across the thernocline vfhich 
is ono-way. The fluid particles ;/hich wove upvfard from the region beloi; 
the thernocline are rapidly deforned and laixed tiiroughout the upper layers; 
fluid particles ;.'hich nove do;:mm'd tovfard the quiet layer below the KLD 
are buoyed upward intact. Thus, the upper layer increases in thickness, 
but decreases in temperature, it the expense of the lo\;er layer. This sane 
phenonenon is observed when the turbulence is induced by v/ave action or in- 
ternal wave action [l] . 
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